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TECHNICAL MEMORANDUM 1417 


ON THE SPECTRUM OF NATURAL OSCILLATIONS OF 
TWO-DIMENSIONAL LAMINAR FLOWS* 

Ey D. Grohne 

1. INTRODUCTION AND STATEMENT OF THE PROBLEM 


In the investigation of stability of a two-dimensional laminar flow 
with respect to small disturbances, we describe a disturbance of the 
stream function moving downstream (in the direction of the x-axis) by the 
"partial wave formula" 


* = cp(y)e ia ( x - ct ) (1.1) 


and obtain then for the distribution of the disturbance amplitude cp(y) 
at right angles to the main flow the so-called stability differential 
equation of the fourth order 

(U - c)(q>" - <x 2 cp). - U"cp = - 2a 2 <p" + <x%) (1.2) 


where U(y) designates the velocity profile of the basic laminar flow 
3h addition, we enforce certain boundary conditions, in the specific 
case of the parallel channel 


cp(±l) = 0 cp* (±1) = 0 


(1.3) 


*"Uber das Spektrum bei Eigenschwingungen ebener Laminar stromungen." 
Zeitschrift fur angewandte Mathematik und Mechanik, vol. 34, no. 8-9, 
August-September 1954, PP* 344-357- 
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which express the fact that even the disturbed flow adheres to the 
bounding walls. In these equations, the velocities and c are 
referred to a velocity of reference furthermore, the lengths x, 

y, and l/a to half the channel width' b, ancl finally jthe time t to", 
the time unit b/UQ. The Reynolds number R r is defined by 


R 


2£ 

V 


The boundary-value problem consisting of differential equation and 
boundary conditions determines, for each pair ^ of parameters a and- R, 
a spectrum of an infinite number of eigenvalues c n . The associated 
disturbances (l.l) are damped when Im( b n ) < 0, and are excited when 


Im(c n ) > 0; a is assumed to be positive apd_real. A basic flow is 
called stable for a value of R when the entire eigenvalue spectrum 
c n , for all possible values of a, contains only dampe d^Hi sturbances . 
Thus the range of the Reynolds number $ is divided up "Into a region 
of stability 0 < R < R* and a region cj>f instability R^> R*, which are 
separated from one another by the stability boundary R*. 


Since, in the literature published .up till now almost exclusively 
neutral oscillations - at most, excited .oscillations- - have been investi- 
gated, we shall investigate in the present report, following a suggestion 
of Prof. Dr. W. Tollmien, the entire spejctrum qf the eigenvalues e n as 
a function of a and R; for simplification, we shall emphasize the 
dependence on aR. A general solution ojf this problem is possible in 
the following two special cases: (l) in; the case U = const, which is 

equivalent to U = 0. We deal here with the "oscillations of a fluid 
at rest" already treated by Lord Rayleigh. The solution is possible in 
the domain of elementary and transcendental functions. The second 
special case concerns the rectilinear Copette flow U = y_ Investigated 
by L. Hopf (ref. 5)- The solution can be reduced to tabulated Bessel 
functions . ’ 


For more general velocity profiles :U(y), the eigenvalues c n can 
be determined approximatively analytically in the following limiting 
cases : 


I. 

In the limiting 
e j.genvalues < 

case 

: n 

aR — > 0 

for arbitrary order n of the 

11 . 

In the limiting 

case 

n oo 

for constant aR 

III. 

In the limiting 

case 

aR — » oo 

for resin 1 ic ted order n. 
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A continuous transformation of the three cases into one another 
for constant subscripts is possible in the above named special cases 
U = 0 and U = y. The assignment of subscripts of the eigenvalues 
can be made in the cases I and II according to increasing damping, that 
is, according to the rule ■ — - 

^( c n+l) = ^ ( c n) U-M 

However, this rule is not always applicable to the case III when the sub- 
scripts used are to remain constant for continuous variation of a and 
E. 


The boundary- value problem formulated in (1.2), (1.3) is, generally, 
not self-adjoint; thus, the reduction to the well-known statements and 
estimates of the Sturm-Liouville theory is eliminated. The eigenfunc- 
tions generally do not form an orthogonal system. They do form, however, 
as 0. Haupt (ref. 3) has shown, under certain assumptions, a system of 
functions that is complete with respect to each of the functions which 
satisfy the boundary conditions (1.3) and are four times continuously 
differentiable. This system of f unctions can be transformed into an 
orthogonal system. 


2. THE LIMITING CASES aR -» 0 FOR ARBITRARY ORDER n OF 
THE EIGENVALUES c n , AND a -m» FOR LIMITED dR 


As already found by Lord Rayleigh (ref. 8), the entire system of 
eigenfunctions and eigenvalues in the case of the basic flow U = 0, 
that is, for a medium at rest, can be given as a closed system. Since 
these eigenvalues are suitable for approximative representations in 
the case of more general basic flows also, we shall derive them here 
briefly. In the case of the basic flow U = 0, the stability differential" 
equation (1.2) is simplified to — — 


- 2a 2 4" + a* 1 ® + iaR • c(<&" - a%) = 0 (2.1) 

where we shall denote the eigenvalues by C, to distinguish them from 
the eigenvalues c of the general stability differential equation. 

The equation is solved by each of the functions -- 


$l(y) = cosh o>y - cosh ay ^ ®jj(y) = sinh coy - sinh ay 


f 2.2) 



4 


MCA OM 1417 


if we put - - 

03 2 = a 2 - iaR • C ' (2.3) 

The part $(±l) =0 of the boundary conditions* is identically satisfied. 
The remining boundary conditions 4>'(±l) = 0 lead to the related branches 
of the eigenvalue equation: -- ~ 


a tanh a = tanh o> in the case I 

\ (2.4) 

a coth a = oa coth od in the case II J 

The equations (2.4) have, for positive q, no fpots o> outside the 
imaginary axis of the complex co-plane. With a>, also -co, is an eigen- 
value associated with the same eigenfunction . Thus it is sufficient 
to consider only the positive imaginary Eigenvalues oo. „If we designate^ 
even eigenfunctions by even subscripts arid odd 'eigenfunctions by odd 
subscripts, the equation (2.4) may have the solutions cdq, o^, cd|j_, . . . 

in the case I and co^, oo^, ... in the case II. The eigenvalues 

can always be made to form a monotonic sequence ' _ 


0 <^<^<^< . : . 

The associated eigenvalues C n are according to equation (2.3) 


( 2 . 

They are, therefore, arranged in the order of increasing damping. The 
nth eigenvalue my be estimated upward arid downward by 


a 2 + * 2 (g + gj 2 < iaR . C n < b, 2 + rt^g + l) 2 " (2.6) 

From the representation (2.5), it follows that the eigenvalues C 
become very large for ctE 0 as well as for. n The same behavior 

occurs, also, for more general velocity profiles U(y) because the main 
parts of the stability differential equation (1.2) are then represented 
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by the equation (2.1). We shall now express this train of thought more 

accurately by subjecting the difference c - C of the eigenvalues c 

to a more accurate estimate compared to the eigenvalue C of equa- 
tion (2.1), for more general profiles. With introduction of the differ- 
ential operators 


L[q>] = u(<p" - cAp) - U"9 M[<p] = 9" - a, 2 <p n[<p] = 9^ - 2a 2 cp" + aSp 

( 2 . 7 ) 

the stability differential equation (1.2) may be written in the form 


i 

Correspondingly, equation (2.l) reads 



( 2 . 8 ) 


Utilizing the fact that the operators M and N are self-adjoint, we 
obtain from these two equations the relationship 


(c - c) J Cfw[$]dy = J $l[cp] dy (2.9) 

If the normalization which is still open for cp is fixed by the rule 



( 2 . 10 ) 


there follows, after introduction of the auxiliary quantities 


Q = 




( 2 . 11 ) 
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from equation ( 2 . 9 ) the representation li. . . . — 

c - c = q + q. ' ( 2 . 12 ) 

In this equation, C and Q may he regarded as known"by virtue of the 
functions $ represented in equations; (2.2). The eigenvalues C n 

have already been delimited in expression (2.6). For we obtain 
directly the estimate 

Q - | f + ^ U • dy + | co | » 1 (2.13) 

In connection with a simultaneous estimate of the function 

(cp - <&)" - <x^(cp - $), we obtain for q‘ ' ■■ - ■ - 

q = for SB ^0 - (2.14) 

\ CO / CO 

If we substitute both into equation (2.12), we obtain, “With consideration 
of equation (2.5) aud expression (2.6), the two partial statements • 


‘n - c - ■ t - °(rTl) 

=n - C n > | f* U • dy + o(f ) 


for acH — > 0 and for arbitrary order n 

1 j2.15a) 

.IT" TT 

for n -> °° Tor fixed JaR 


(2.15b) 


The latter estimate- indicates that the ‘eigenvalues : of the stability 

differential equation for sufficiently high dfder n finally tend toward 
the eigenvalues C n of the "zero flow" (with the real ’part- increased by 

the mean velocity of the basic flow). (Compare eq. (2.1).) A mutual 
coordination of the eigenvalues c n to i the eigenvalues i_£J n , however, 

is by virtue of equation (2.15b), meaningful "only due ta.the fact that 
the difference j C n - C n+l! of the approximation eigenvalues comes out 

considerably larger than the estimated remainder in equation (2.15b): 

For, because of (2.5) and (2.6) 
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c n " °n+l = constant * (2.16) 

is valid. It should he mentioned that F. Noether (ref. r J, p. 239 , for- 
mula (28)) has already indicated an asymptotic representation for slightly 
differently defined eigenvalues for unlimitedly increasing order, 
although only intimating an argument - which leads one to expect con- 
siderable diff iculties . . — . 

We mention, furthermore, an estimate for the eigenvalues c indi- 
cated by C. S. Morawetz (ref. 6, p. 580) 

c - c^ < A • (aR) -1 / 2 

where e n is an approximative eigenvalue which (in our notation) is 
determined by the equation 

m j* ki? - o • * - n 

and corresponds more or less to our approximate eigenvalue C n intro- 
duced in equation ( 2 . 1 ). In the above estimate of Morawetz, neither aR 
nor n may become arbitrarily large; in the first case, the eigenvalues 
would shift into the excluded neighborhood of c = w^y^^w = U; y^ des- 
ignates the wall); in the other case the estimate would become meaning- 
less since the behavior of the quantity A for unlimitedly increasing 
n is not given. 

3 . RECTILINEAR COUETTE FLOW. THE LIMITING CASE oR 

FOR FINITE ORDER n 

In the special case of the basic flow; at rest, U = 0 , the behavior 
of the eigenvalues for unlimitedly increasing aR is described 

by the formula (2.5) in which the quantities cc^ no longer depend on 

aR. In deviation from this law, there results for more general velocity 
profiles a behavior like ; " " ; 
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r . 

C n - u(-l) = -=-£- for . aR -4 00 (3.1) 

Xjai - -fr 

where the complex valued quantities no longer depend on aR. If _ 

these eigenvalues are adjoined to the : eigenvalues in equation (2.15a), 
by continuous transition of a and aR, the ordering principle (1.4) 
according to increasing damping is lost even in special cases like the 
rectilinear Couette flow. If we therefore desire that, the subscripts 
of the eigenvalues remain unchanged for continuous transformation 

of the limiting cases aR -» 0 and aR -»» into one another, we must 
actually carry out this procedure which presupposes a^general solution 
of the eigenvalue problem or a solution which is approximate only inso- 
far as the individual eigenvalues still remain distinguishable from one 
another. We succeeded in obtaining a! solution in this^ sense only in the 
special case of the rectilinear Couette flow. It will, therefor e, form 
the subject of the following section. 

After insertion of the velocity profile U = y of the rectilinear 
Couette flow into the stability differential” equation _(l. 2), the latter 
can be reduced, by means of the substitution 

♦ = 9" - a 2 ? _ (3*2) 

to the Bessel differential equation in the auxiliary form 

\|r" - [ictR(y - c) + a^J \|r = 0 (5-3) 

In order to arrive, through the boundary conditions, at the eigenvalues, 
we must invert equation (3*2) in the form “ 

*y)- f* *(i) sl 4 «/ . - 1) an 7" (3-M 

The boundary conditions <p(-l) = cp'(-l) -0 ^then are identically satis- 
fied; the remaining boundary conditions cp(+l) = 9'(+l) = 0 require 
that the two equations . 4: -rr • 



sinh ay 
a 




^(y)cosh 


ay dy' = 0 


(3-5) 


hold. 
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* 

By means of the substitution 


y - y, = 6 T] with e = (aR)" 1 /^ and y. = c + ( 5 * 6 ) 

*■ *■ coR 

the differential equation ( 3 * 5 ) niay be transformed into the differential 
equation 


If 


i + Vtf - 0 (3-7) 

*7 


i|f Tj) and * 11 ( 1 ) 


(3-8) 


are assumed to be two suitable fundamental solutions of this equation, 
and T)_-^, is assumed to designate the values which are, because of 

equations ( 3 . 6 ), associated with the walls y = -1 and y = +1 


'Ll = I l 1 " c - 


±af 


■K 1 - 


±QL C 


aR / ^1 € V" ~ aR 

there follovs from equations ( 3 - 5 ) the eigenvalue equation 


(3-9) 


+ x C'n) - lD | r € - T) d 7 * J^1_ 1 + n (T i )cOBh aeT i dT > - 


V o\l 


f ^n( 1 ) S1D ^ £ a€ " dT l • F * (T))cosh aeTj dr\ = 0 (3*10 ) 

j Ti T]_-^ 


7-1 


For further treatment of this equation, the introduction of a 
sequence of functions q) by the Laplace integral 


An( 7) 


=— r 

2rti J a 


r|Z+i 


z 11 - 1 . dz 


(3.11) 
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is advisable in which the path of integration" A runs from infinity 
to infinity in the manner drawn in figure 1. . The functions A n ( q) 
satisfy the differential formula 


dr, 


= Am-iU) 


and the recursion formula 


(3.12) 


i • + rj • A n+1 + n • Ajj^ = 0 (3*13) 

by means of which all the functions A n Cr() and their integrals and 
derivatives may be constructed recursively from the three basic functions 


Aq( t)) a-|_( rj) A^Cn) (3.1 1 *-) 

The significance of these functions: A^t]) for the stability prob- 

lem lies in the fact that the two particular solutions of. the differential 
equation (3-7) needed in the eigenvalue equation (3.10jcan be repre- 
sented in the form r " 7 


2*i • - / gjri \ 

t-j-Cn) = A 1 ( q) ^ I3 -(t]) = e : 5 • A 1 U . e 5 . ) (3-13) 

That the differential equation (3*7) is satisfied follows from the 
formulas (3.12) and (3.13). The linear independence of the two functions 
follows from the fact that the Wronski determinant, which is constant of 
course, does not disappear at the point T) = (TT 


The basic functions A^(r|) and Agft]) = A^'(t)) have been numeri- 
cally tabulated (in somewhat different nbtatioh) for a quadratic point 
grid with the mesh width 0.1 within the bircle 1 I tj | = 6 _ of the complex 


tj-plane by H. H. Aiken (ref. 1). 


The basic function 


A°(n) - J A^TjJdr] 


can be determined from it by a numerical integration. 


Outside of this table , the behavior .of the functions A^rj) may be 
inferred from the asymptotic series representation 


JM. *■*** -f * 

J nv 


v=0 


(3-16) 
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which is valid for 1 7 | -> 00 in "the angle space 

with arbitrarily small 8 > 0 . According to H. 
first coefficient of the series is 


- ^ + 8 = arg 7 = - 

O 6 

Holstein (ref. 4 ), the 


* 11,0 ~ 




(3.16a) 


The asymptotic series are obtainable directly from the Laplace integral 
(3*11) by means of Riemann's saddle-point method. 

For the representation of the eigenvalues, the zeros 7^ of the 
function Aq(t|) are necessary. An asymptotic calculation of these zeros 
for |%| » 1 is no "k possible directly by means of expression (3*l6), 

since the zeros would move out of the range of validity of this repre- 
sentation. VTe avoid this difficulty by applying the second relations 
obtainable from the integral (3*11) 


-i 


Jt 


iitn 


AnVl e Z ) - e 3 


A*Vf * - 

(* = con jugate- complex values) 


i*\ 

T 


( 3 - 17 ) 


and 


i £jm / ± 2rt\ ± 4rtn / ± 4*\ 

A n (7 1 )+e 5 A n (^e 3 ) + e 3 A^rje 3 / = P n (r 1 ) ( 3 . 18 ) 

where the polynomials P Q ( 7) of the degree -n satisfy the same recur- 
sion formula (3*13) ^ 


■*- ‘ p n+3 + lPn+1 + ^n ” ® 
with the initial elements 


( 3 . 19 ) 


P 0 = 1 


?1 = 0 P 2 = 0 


Combination of formulas (3.16) and (3.18) then yields the asymptotic 
representation valid for I z I in the angle space I arg z) < n - 8 
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f i co Jz z 3/2 + *\ 

AoW & )il- :V. ■■ /k k ' (5.20) 

}Jj( . Z^' 

Hence, there follows, for the zeros of A^q^) = 0 which, according 

to equation ( 3 • 17 ) 5 lie in pairs symmetrically with res pe ct to the 
straight line arg = 5 jt/6, the asymptotic representation' 


% = (3*n) 


2/3 


i 

e *- 


5* 4. 1n(n»2^\f3 ) 
3im ■ J 


n = 1, 2, 3 

The value of the lowest pair of zeros was calculated to he 


( 3 - 21 ) 


% - ^-257 


iPp ±0.2708 

e L -= - 


according to. the table of Aihen (ref. l)i. 


(3.21a) 


For further treatment of the eigenvalue equation (3_.10), it is 
advisable to expand the functions sinh qeq and cosh ccTeq into their 
Taylor series, and then to interchange the summation with^ integration 
which is justified by a theorem Of Bromwich. (Compare ref. 2, p. 398.) 
The series obtained converge, according to theorems of the Laplace trans- 
formation, for each value of ae. If these series are broken off after 
the first terms, provided with residual terms, 'and substituted into the 
eigenvalue equation (3*10), the latter is, forhhis reason, and with 
consideration of equation ( 3 • 17 ) > simplified to - 


M’i) ' A o( ,, -i) 
M\) - A 2C-l) 


*SK) - ^(-<1) 
A IK) - A l(-fi) 



for 



= constant 


(3.24) 


What happens now when aE increases; beyond all limits, that is, 
when e tends toward 0? Because of the. relationship following from 
equation ( 3 . 9)5 \ ~ ^ 1 = at least; one of the two quantities 

q_^, q ^ must tend toward infinity for k -*0,1 on a parallel to the 

real axis. It is sufficient to assume, this regarding q^, because in 

the other case everything would form a mirror image with respect to the 
imaginary axis of the:, q-p lane. (as essentially occurs in the Couette flow 
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where with C, also, -C* is an eigenvalue). With consideration of 
the asymptotic "behavior (eq. ( 3 • l 6 ) ) , the eigenvalue equation (3.24) 
then is simplified to 


i — - £ = ofe 2 ) for e « 1 

^(l-l) 2 . 


(3-25) 


From this formula, we recognize that rj for e 0 must tend 
toward the zeros of the function A^( t]) estimated in expres- 

sions (3*21). We thus obtain for the eigenvalues c, with considera- 
tion of equations (3-9)) the asymptotic representation 


c + 1 = - € % + o(e 2 ) 

( 3 . 26 ) 

with from k Q ^ = 0 

Thus we have proved the previously given eigenvalue formula (3*1) 
for the special case of the Couette flow. 

In order to follow the variation of the eigenvalues c over the 
entire range 0 % aR < 00 , we must go back to the eigenvalue equa- 
tion ( 3 - 10 ) or its approximate form in equation (3-24), with the func- 
tions Aq(tj), A^(q), A^( t]) to be assumed as known. We have accord- 

ingly calculated the 12 lowest eigenvalues c as functions of oft for 
a fixed value a = 1 and represented them in figure 2. The variability 
of the eigenvalue curves with a is only slight and becomes, for 
instance, for aR -» <» with e small of the order O(a^e^). 

The eigenvalue spectrum of the rectilinear Couette flow has been 
discussed already by L. Hopf (ref. 5). Hopf replaced, more or less on 
the level of our eigenvalue equation ( 3 - 10 ), the solutions \|r^, i|r__ 

represented by him by Hankel functions of the order 1/3 - by the first 
terms of their asymptotic series ( 3 . 16 ), whereby the eigenvalue equa- 
tion was simplified to an algebraic equation of auxiliary arguments and 
circular and hyperbolic functions. However, since Hopf committed certain 
errors in the asymptotic representations of the Hankel .functions, his 
results require partial corrections. Although these changes are hardly 
significant for small values of aR, the values of, for instance, 

in a formula corresponding to ( 5 . 2 6 ) undergo a considerable change. The... 
topological connection of the eigenvalue curves c = c(a,aR) also 
appears different to us from what it appeared to Hopf. However, the' 
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qualitative picture of the eigenfunctions , the physical^ conclusions 
drawn from it, and the main result - thpt all oscillations are damped 
remain the same. 


4. THE LIMITING CASE oR -> » FDR FINITE ORDER "n FOR 
SYMMETRICAL BASIC FLOWS 


For a basic flow with symmetrical velocity profile 

U(y) = U(-y) (4.1) 

the stability differential equation (l.;2) always has a fundamental 
system of four solutions cpy <f> 4 so that 


and 


9 2 (y) 

\(y) 

are even functions of .y 

\(y) 

9,(y) 

3 

are odd functions of y_ 


-(4.2) 


If a linear combination of these solutions is_ to satisfy* the boundary 
conditions of-equations (l.3)'in the sdquence_ cp( -l) = 0, 9 f (-l) = 0, 

cp(+l) = 0, cp l (+l) = 0, the . following determinant, simplified with 
consideration of the symmetries (4.2), must disappear: 



^(-D 

$ 3 (-l) 

?£(-!> . 

fy(-l) 

?{(-!) 

^(-D 

^(-1) .. 


-^(-D 

-9 5 (-l) 

^(-1) , 

9 4 (-l) 

${(-1) 

^(-1) 


-?i(-l) 


Since this determinant may be written &s therproduct of. the two-column 
deterrdnants - . 


9 X (-1) 

? 3 (-l) 


<p 2 (-i) 

? 4 (-i) 

9;(-i) 

3 


^(-i) 

*;<-*> 


(^• 3 ) 
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one obta ins , by equating one of the two factors to zero, one branch of 
the eigenvalue equation each time. For this reason the eigenfunctions 
can be either only even or only odd, with the respective eigenvalues c 
generally being different. 


In order to arrive from these equations at asymptotic eigenvalue 
formulas, we shall determine the four fundamental solutions (4,2) 

in such a manner that they are available for appropriate 

asymptotic expansions. We find that the fundamental solutions described 
by W. Tollmien (ref. 12), "Asymptotic Integration of the Stability Differ- 
ential Equation", the asymptotic representations of which are provided 
with residual- term estimates, are suited to the problem. 

In order to establish the connection of these fundamental solutions 
with ours, it is indispensable to discuss first the concept of "friction- 
less approximation." The quest for solutions of the complete stability 
differential equation (1.2) which for dR «, together with their 
derivatives with respect to y, tend toward a limiting function 

lim cp(y,dR) = x(y) (4.4) 

aR-»oo 

leads to the so-called "frictionless differential equation" 

(u - c)(x" - a 2 x) - U"X = 0 (4.5) 

which mast necessarily be satisfied by such limiting functions. If we 
want to use the solutions of this frictionless differential equation 
for the approximation of the solutions of the complete differential 
equation for aR ->», we must not disregard the range of validity of the 
boundary-value statements in equation (4.4) in the complex y-plane. 
According to W. Wasow (ref. 13), the following theorem is valid with 
respect to this; 

"Of the four fundamental solutions (4.2), one even and one odd 
solution can be determined in each case so that with two appropriate 
frictionless solutions 3q_(y) and Xg(y) the approximations 


h (y) = XjCy) + o(^) 

<P 2 (y) = yy) + 0 ( 4 ) 


X-^(y) = odd function of y 
Xg(y) = even function of y 


(4.6) 
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in each, fixed interior of a double region (l_+ II) or (II + III) or 
(ill + I) are valid and become invalid) in each case, in the comple- 
mentary third region III or I or II. The same is true jf or the deriva- 
tives with respect to y." (Compare fig. 3«) 

The boundaries between. the regions I, II). and III satisfy the 
equation • - — • — - 



if y denotes the "critical point" defined, by 
4 

u (y k ) = c Re y k < 0 _ (4.7) 

For more details regarding the regions -I, II, “and III see Wasow (ref. 13). 

The frictionless differential equation (4.5) has at the critical 
point, U^y^ = c, a singular point with regard to determinateness. Two 

fundamental solutions take the form ; ~- 

x x (y) = (y - y k ) • P x (y - y k ) . % 

K 

x 2 (y) = R(y - y k ) + or ^(y - y k ) • (y - y k )in(y - y k ) (^*9) 


denote power series with the.. beginning- 


p^z) = 1 + z • + o(z 2 ) P 2 (z) = 1 + o(z 2 ) (4.9a) 


(Compare W. Tollmien, ref. 12, p. 35*) The common radius of convergence, 
of these power series is limited either by the radius of convergence of _ 
a corresponding series for U - c or by the next -adjacent zero of U - c 
as a singular point of the differential equation. (7 


if P and P 
“1 “2 
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For the further development it is advisable to introduce a sequence 
of functions B^Cq) by the Laplace integral 


^ 'ST 

which is comparable to equation (3-11). In it, c = 0.577 2 . . . denotes 
the Euler constant. The path of integration ]3 runs, in the mann er 
indicated in figure 4, in the complex z-plane cut open along ( 0, -i°=) 
from infinity to infinity. 

The functions T) ) satisfy the differential formula 




T]Z+i -j 


z n_1 ( In z + c)dz 


( 1 . 10 ) 


«3B n (n) 

dr) 


= Bn+lC^) 


(l.ll) 


and the recursion formula 


1 ' -^+3 + ^ * ^+1 + n * ®n = p n ( 1 . 12 ) 

in which P n ( rj) are the polynomials defined in equation ( 3 * 19 )* 

By means of these two formulas all the functions B n (Tj) and their 

derivatives and integrals can be constructed recursively from the three 
basic functions 


B 0 (ri) B-lCtj) B 2 (t!) 


By means of the representation 


B n U) = 2*i 

(n=l,2) 


AqU) 


A*(-t]*) 

(-i) n A*(-n*) 


(1.13) 


( 1 . 13 a) 


(* = conjugate-complex value), the basic functions and B^ can be 

reduced to the functions (Compare W. Tollmien, ref. 10, p. 27.) 

The significance of the functions B n (T)) for our stability - eigen- 
value problem lies in the fact that the function B-^(t)), because of 
equations (4.11) and (4.12), satisfies the differential equation 
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d^B, d% 

1 “IT + r » “V = x (^) 

drp dq 2 

which, with the designation "differential equation for _the friction 
correction, " has been introduced as an essential constituent into the 
asymptotic integration of the stability differential equation by 
W. Tollmien (refs. 10 and 12). 

After these preparations, we turn to thefour fundamental solutions 
cpjj^ cpjy of the complete stability differential equation 

constructed by Tollmien, regarding its (ability to be expanded asymp- 
totically. According to W. Tollmien (ref. 12,. p. 77) these four solu- 
tions may be determined, with use of the substitution 


= c (4.i6) 

in such a manner that they have in a fixed interior of the rj-plane (com- 
plex for reasons of analytic continuation) as well as iif every fixed 
interior of the region II of the frictionless_approximation (compare 
eq. (4.6a)) the following asymptotic representations: 

9j(y) = X x (y) + o( e 5) (4.17a) 


'pjjCy) = £2^) + • e i( £t i) * € |b tU) +* i+ n 

k L- 1 ' 

or in every fixed interior of II ' 1 

c p II (y) = *2^ + o( e5 ) 

Furthermore, 

9m(y) = A_ 1 (rj) + 0(e) in.. |t)| = constant - (4.17c) 

is valid. Finally, there applies, according to. W. Waspw. (ref. 13), 
quotient-asymptotically in every fixed interior" of H (compare eqs. (4.6)) 
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yy), yy) ~ constant 



exp 


jiaR 





Corresponding formulas are valid for the derivatives. 


For farther treatment of the eigenvalue equation (4.3); we must 
express the fundamental solutions 9^ ... 9^ used in it by the above 

fundamental solutions 9j • • • 9jy If for the latter, the representa- 
tions (4.17) are used immediately, and with the residual terms in each 
circle jr]| = constant for e -> 0 being valid, the result reads 

+ o(e 2 ln e) 

* 


(4.21a) 

from the frictionless solutions and Xg- (Compare eqs. (4.9)) 

Furthermore , 


U k 


ypCy) "^1,2^111^) “ 1+e ' TT 7 

u k 


i]^ln e+ 1 + + ®-i( *)) 


ujl Xp(0) 

” lth m S 1 - - y 


U£ _ X£(0) 

w s 2 ~ ~ x-Toy 


<P 5 yy) = ^mCy) + °U 


5/4 


constant^ 

“eVF _ 


= Ay i) + 0(e) 


(4.21h) 


is valid. Corresponding formulas are valid for the derivatives. 


If we 
tion (4.3) 


now write the two eigenvalue equations obtained in equa- 
as a product of three factors, for instance. 





?}(-!) 5{( -1) - *P^(-1) 
fjll) " qy-i) - 4 3 (-d 
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( i3 = arbitrary constant), the zeros of .the two first factors do not 
make a contribution to the eigenvalue , configuration since they are 
compensated, £y corresponding poles of the third factor - unless the 
derivative cp^(-l) should disappear simultaneously. It is therefore 

sufficient to find only the zeros of the third factor. After insertion 
of the approximations (4.21) we thus obtain 



(1 - 1_x) 


or 


a q( n) 

a_ x ( n) 


+ eKn) 


2 * e + ^ + ®1,2 + Bo(,)] 

U k .. - . • 

ytt " “ 

1 + e rjT e + 1 + s l,2) + B-lUj 

u k L ' — 


+ 0 


(e 2 ln e) 


(4.22) 

with T)=rj =(-1- y k )/e and S^, according to equations (4.21a) 

from the frictionless solutions. The function \|r( T}) stemming from a 
next-higher approximation in equation (4,21b) reads 


M *^$ii 


A _ 2 (n) + | A-jU) - 

A_ 1 ( Tj) 


10 a 4^ ^ 


(4.22a) 


and may be reduced, by means of the formulas (5. 12) and (3.13)> to the 
three tabulated basic functions Ay, A^, A^. 

How do the eigenvalues c behave !if inTthe * eigenvalue equa- 
tion (4.22) we let aR -» <*>, that is e -4 0? Evidently rj ^ then tends 

toward the zeros t)^ of the function Aq( r]) . By Taylor series expanded 

about these zeros, there follows more exactly 







^(hn) 

■ 


e In e + 0(e) 
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The eigenvalues then "behave asymptotically like 


c _ u(-l) = -U£ . €T! n 


U k 

U k 


, A 2 (%") 
A l{%) 


e 2 ln e + 



(4.23) 


with 1 ^. from A^ =0. As a supplement to equations (3.21) we shall 
give here a few zeros t}^ and values iAg/A^: 


% 


“2 (%)/%(%) 


-4.122 + i . 1.065 

- 1.686 - 

i . 1.222 

- 2.983 + i . 3.037 

+I .902 + 

i . 0.851 

- 6.8 + i . 2.5 

- 2.2 

i . 1.5 

-5-5 + i . 4.5 

+2.4 + 

i . I.l4 


For | % |»1 
there applies 


M%)_ 

A i (%) 


(4.24) 


The remarkable fact about the asymptotic eigenvalue formula (4.23) 
i 6 that it is transformed into the corresponding formula ( 3 * 26 ) after 
substitution of the velocity profile U = y of the Couette flow, although 
the two formulas were derived under completely different assumptions. 


The asymptotic eigenvalue formula (4.23) is already so greatly 
reduced that it no longer permits a distinction of the eigenvalues c 
which are associated with even or odd eigen functions. For this, we 
must go back to the more exact formula (4.22) in which the character- 
istics "even" or "odd" of the eigenfunctions are taken into consider- 
ation by means of the constants S-, and S Q , to be determined "without 
friction." 2 


We have used the eigenvalue equation in the form (4.22) also for 
the numerical calculation of the eigenvalues c in the examples treated. 
We selected as examples the two-dimensional Poiseuille flow and a flow 
with an inflection-point profile. We represented the variation of the 
four lowest eigenvalues as functions of R for a fixed value of a 
in figures 5 and 6 . The numerical calculation itself is - after reduc- 
tion of the nonalgebraic elements contained in the eigenvalue equation 
to the three tabulated basic functions A^, A-., A^ and to the fric- 

tionless solutions - a problem involving numerical methods, the details 
of which ca n not be discussed here. We shall mention only the following 
approximate representation of the frictionless constant ■ 
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— . So - A . a" 2 + 0(1) for a <£< 1 with A 
Ui 2 ..... 

k 


(4.25) 

(Compare W. Tollmien, ref. 11, p. loo), which may he applied advanta- . 
geously for small values of a. 

The subscripts for the eigenvalues c"" obtainedlTrom equation (4.22), 
in the sense of a continuous connection witfi the limiting case dR -> 0, 
remain an open problem here. In the range qf validity of equation (4.22) 
alone, a generally valid choice of subscripts according to the rule 
Tm (c^i) Im^c n ^, that is, according to increasing damping, cannot in 

principle be carried out, either. The! zeroi^in equation (4.23) can be 
ordered according to the increasing imaginary part, but the Im(c n ) 
curves may penetrate one another if e is changed. ~ ^ ' 

5- THE FRICTIONLESS EIGENVALUES WITHIN THE LIMITING CASE dR -> » 



J -1 


Determination of the Excited Eigenvalues 

Let the approximation (4.6) by means, of _the frictipnless solutions 
be suited either to the double region % + 11^ (compare fig. 3) or to the 
double region II +-III whereby the logarithmic term is .always uniquely 
determined in the frictionless solutions. Applying the approximations 
(4.17), we then obtain, by way of the eigenvalue equations (4.3), eigen- 
values c which, for dR °°, tend toward the so-called "frictionless 

eigen values" c^^a) which are defined by^he boundary condition 
-1) =0 or Xg( -1) =0 of the odd :or even frictionless solutions 
Xp, Xq. The following general statements may be made regarding these 

frictionless eigenvalues, limited by th;e r.ange of validity of the 
boundary-value expressions (4.6), according tQ W. Tollmien (ref. ll), 
partly on the basis of the "Ray leigh-Tollmierr- theorems 

"For velocity profiles without turning points,- no excited friction- 
less eigenvalues are possible. The approximation. (4. 6X associated 
with the damped frictionless eigenvalue's ^must always take place in the 
interior of the double region i + II." 

"For inflection-point profiles, there always exist excited friction- 
less eigenvalues associated with an evert eigenfunction." 
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Beyond these general statements, frictionless eigenvalues associ- 
ated with an odd eigenfunction were not found in any of the examples; 
neither did we find eigenvalues such that the associated approxima- 
tion (4.6) would have taken place in the interior of the double region 
II + III. As examples, we chose the two-dimensional Poiseuille flow as 
representative of a profile without an inflection point, and the inflection 

point profile U = ^\[2 - l) + ^2 - x cos y. The frictionless 

eigenvalues c, found only associated with an even eigenfunction, are 
represented in figures 7 and 8. The range of existence of these eigen- 
values is always given by an interval 0 = a = constant. For the fric- 
tionless eigenvalue c, Tollmien (ref. 11, p. 100) has set up the 
following approximate formulas : . _ 





with U(y K ) » c r 
(5-1) 


We now seek the connection between the frictionless eigenvalues 
and those discussed up till now. The closed solution, in the case of 
the Couette flow, cannot give an answer to this problem because the 
frictionless eigenvalues in question do not exist there at all. How- 
ever, it is possible to insert the frictionless eigenvalues into the 
equation (4.22) and thus to interpret them as a limiting case within 
the eigenvalues (eq. (4.23)). 

Let us, therefore, perform on the eigenvalue equation (4.22) the 
limiting process aR °°, that is, e -> 0 for constant ( -1 - y^.) = eri ^ 

the justification of this procedure is based on the equation preceding 
(4.22). By means of the asymptotic formulas (3-16) there then follows 
for aR -» w 


1 + (- 1 - y k ) ut(M-i ~ yk) + s i,2*l “ 
k 4 J 


These are, however, precisely the first Taylor terms of the frictionless 
eigenvalue, equation X]_(-l) =0 or Xg(-l) = 0 which would be obtained, 

according to the significance (equation (4.21a)) of -3^, Sg, in the case 
of Taylor expansion in the sense of the series (4.9). 
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On the "basis of this findings the determinationjpf the excited 
eigenvalues (imc > 0) can be simplified. Since, according to the results 
of the second section, excited eigenvalues can appear only within the 
first eigenvalues of finite number ahd for sufficiently large values 
of R, it suffices to examine equation ( 4._22 ) with resp ect to excited 
eigenvalues. For' this, the following alternative is _valid: Excited 

eigenvalues can be (approximately) determined either .by the friction 7 
less boundary-value problem in combination^with sufficiently large values 
of R, or they lie in a neighborhood’ of cT= U(-l) and can be determined 
by means of one of the equations (4.22) or (4.23) as associated with 
finite values of t] 

As follows from this for sufficiently large values of R, but as 
was confirmed in the examples for the smaller values .of R also, the. 
greatest excitation for inf lection-point profiles is .always combined with 
the frictionless eigenvalue or its continuation toward smaller R 
values. Hence, there follows the well-known fact that, in the case of 
turning-point profiles, the stability behavior may be concluded even 
from the frictionless differential equation alone. Let us compare to 
this the calculation of the frictionless eigenvalues of G. Rosenbrook 
(ref. 9) for an inflection-point profile which he had, . measured, in a diver- 
gent channel. ' ■' Ti —■ 


6. FORM OF THE EIGENFUNCTION. THE INNER FRICTION LAYER. 
THE VARIATION OF A DISTURBANCE WITH TIME. 


In order to judge the variation with time of a disturbance, we shall 
decompose the latter into partial waves of ’The type In equation (l.l). 

It is then necessary to know the variation of the amplitude cp(y) over 
the channel width. We consider here only the case of very large values 
of ctR. 

For the Couette flow, there follows, by equation^ 3. 4), in the 
notation of equations ( 3 - 6 ), ( 3 - 9 ); and ( 3 . 13 ) for e « 1 , that is, 
aR » 1, as approximate expression for the eigenfunction 


9(y) = F( l) 



with 

F ( ) . A -i (T>) ; Mvl) 

^ " A-i(Ti) - 


( 6 . 1 ) 


( 6 . 2 ) 
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The bo undar y conditions cp( — l) = 0 are identically satisfied) the 
remaining boundary conditions are identical with the eigenvalue equa- 
tion (5.24) . As follows even from the differential equation alone , 
cp*( -y) is an eigenfunction associated with the eigenvalue -c*. We 
calculated accordingly for a = 1 and R = 10^ the eigenfunction 
associated with the eigenvalue c = - 0.7 - i 0.3 and represented it in 
figure. 9 * 

It is striking in this figure that the essential changes of the 
eigenfunction occur in a layer -1 = y = y^ which could be defined 

perhaps by the angle space arg tj = Jt /6 of the strong increase of • 

A ^(tj). In the variable y this "inner friction layer" is, according 

to equations ( 3 * 6 ) and ( 3 * 2 6 ), approximately 

- 1 ^ y^ -1 + e^-Re + Im e«l ( 6 . 3 ) 

This representation shows that the width of the layer increases with 
growing order n of the eigenfunctions; the magnitude of damping"' ~~ 
increasing simultaneously. The velocity of the associated disturbance 
wave is approximately equal to the velocity of the basic flow in the 
center of the layer. Furthermore, the thickness of the layer tends with 
e toward zero. The physical interpretation of this situation signifies 
according to Eopf (ref. 5 , p. 57 ) "that any arbitrary disturbance for 
large values of R is damped in such a manner that, finally, disturbances 
seem to emanate only from the walls, without mutual interference - a 
behavior which reminds one of frictionless fluids." 

For more general basic flows. Toll mien (ref. 12) set up an approxi- 
mate expression for the eigenfunction; in it, one can recognize again, 
in the case of damping, an "inner friction layer" which would have to 
be defined by the angle space jt /6 % arg tj ^ 5 if /6 of the great changes 
in increase of B_-j_( q) or A_^( rj) . In the variable y, this layer is, 

according to equations (4.l6), 


-1 + 


c r + ^ c i < < 

u’(-i) y 


1 + 


c r “ ^ c i 
U’(-l) 


€ — > 0 


(6.9) 


whence we obtain for the higher eigenvalues, according to equation (4.23), 
again the formula ( 6 . 3) for the Couette flow. If, however, frictionless 
damped eigenvalues c in the sense of section 5 exist, the inner fric- 
tion layer, expression ( 6 . 9 ), retains also for the limiting process 
e ->-0 a finite thickness and a finite distance from the wall. We calcu- 
lated, for this latter case, the even eigenfunction in the example of 
the Poiseuille flow, for a = 1 and R = 7-7 X 10^, and represented it • 
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in figure 10. The eigenvalue c = 0.178.- i X 0.0^9 hardly deviates 
from the frictionless eigenvalue associated with a = 1. 

Comparing the inner friction layer with the boundary layer, we may 
say that the boundary layer represents that flow region in which the 
behavior of the laminar basic flow is decisively influenced by the inner 
friction, whereas the inner friction; layer. . indicates the region where 
the disturbance is decisively subject to the influence of the friction, 
since outside this layer the disturbance can be determined without 
friction. ; 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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Figure 7. - Two-dimensional Poiseuille flow. The frictionless eigem 
value c associated, with an even eigenfunction , as a function 
of a . 


4 



Figure 8. - Inflection-point profile. U = (\J2, - l) + ^2 - \/2)cos y ijjl, 

The frictionless eigenvalue c associated with an even eigen- 
function, as a function of a. 




Figure 9. - Rectilinear Couette flow. Eigenfunction <f> '(y) for 
a a l and R = 10® associated with the eigenvalue 
c = -0.70 * i x 0.30. _ 



Figure 10.- Two-dimensional Poiseuille flow. Eigenfunction <p’(y) 
for a = 1 and R* 7.7 x 10^ associated with the eigenvalue 
c= 0.178 - ix 0.049. : — • - 



